Dimerized and trimerized phases for spin-2 Bosons in a one-dimensional optical lattice 
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We study the phase diagram for spin-2 bosons loaded in a one-dimensional optical lattice. By 
using non-Abelian density matrix renormalization group (DMRG) method we identify three possible 
phases: ferromagnetic, dimerized, and trimerized phases. We sketch the phase boundaries based 
on DMRG. We illustrate two methods for identifying the phases. The first method is based on 
the spin-spin correlation function while in the second method one observes the excitation gap as a 
dimerization or a trimerization superlattice is imposed. The advantage of the second method is that 
it can also be easily implemented in experiments. By using the scattering lengths in the literature 
we estimate that 83 Rb, 23 Na, and 87 Rb be ferromagnetic, dimerized, and trimerized respectively. 

PACS numbers: 67.85.Fg,03.75.Mn,67.85.Hj 



Cold atomic gases have been actively studied in re- 
cent years for they offer us new possibilities of study- 
ing quantum many-body systems In the very early 
experiments of dilute Bose gas in a trap, Bose-Einstein 
condensation was beautifully observed directly 0. We 
have also witnessed the realization of the Bose-Hubbard 
model and the observation of the supcrfluid-Mott transi- 
tion 0], a phenomenon theoretically predicted long ago 
but only observed recently. In this case, the presence 
of a lattice and the interatomic interaction actually de- 
stroy the superfiuid, resulting in an "insulating" state. In 
Q and also the many following experiments, the Bosons 
are spin-polarized and so effectively spinless. However, 
Bose-Einstein condensation of bosons with spin degree 
of freedom has also been realized Q . Hence it is natural 
to ask what would be the spin ordering of such Bosons in 
a lattice in the Mott-insulating state, where even though 
one is confined to integer number of particles per site, the 
spins of the neighboring sites can still interact via virtual 
tunneling. Indeed, this question is of much theoretical 
interest, as it can be easily shown that the effective spin- 
Hamiltonians one can realized for spinor bosons loaded in 
an optical lattice are very different from the Heisenberg- 
like Hamiltonians much studied in electronic systems [5[ . 
Similarly, a two-component Bose system allows us to re- 
alize the XXZ spin-1/2 model Q much discussed in the 
theoretical literature. 

In this paper, we consider spin-2 Bosons in a one- 
dimensional lattice. Spin-2 systems are already available 
and experimentally studied |7H1C| . The theoretical phase 
diagram of spin-2 condensates is a function of the scat- 
tering lengths ag in the spin S = 0, 2, 4 channels (Til [T^| . 
It is divided into three regions, which are named ferro- 
magnetic (F), polar (P) and cyclic (C) in Ref. [llj. For 



spin-2 Bosons with one particle per site in a higher di- 
mensional lattice in the insulating phase, it can easily 
be shown that again the phase diagram is divided into 
three regions in the mean-field limit, in direct analogy 



to the Bose-condcnsed case [13Hla ] (see also [17j). In 
one-dimension (ID), however, strong quantum fluctua- 
tions are expected to substantially modify the phases. In 
particular, the polar and cyclic phases are no longer ex- 
pected to be stable. These states break rotational sym- 
metry, which implies the existence of linear Goldstonc 
modes, and thus have diverging quantum fluctuations in 
ID. 

In this paper, we use non-Abelian density matrix 
renormalization group (DMRG) method to determine the 
general phase diagram for one particle per site. This 
regime is much more stable than multi-particles per site 
and thus has a much better chance of being realized ex- 
perimentally We find three phases, which are ferromag- 
netic (F), dimerized (D) and trimerized (T) phases. The 
ferromagnetic state has a macroscopic spin, and has large 
degeneracies arising from the choice of the spin projec- 
tions. The dimerized phase has spontaneously broken 
lattice symmetry, with unit cell consisting of two lattice 
sites. The ground state is a spin singlet, with finite gaps 
to the first excited states. The trimerized phase is the 
most intriguing. At finite size with total number sites 
N multiples of 3, the system is gapped with a spin sin- 
glet ground state. The gap however approaches zero as 
N — > co, resulting in a gapless phase. In this N — > oo 
limit, there is also no broken lattice symmetry. In the fol- 
lowing we shall discuss this general phase diagram and 
present numerical evidence leading to our claims. We fur- 
ther discuss the expected ground states for some available 
spin-2 elements and how the dimerized and trimerized 
phases can be obtained experimentally as well as being 
tested. 

We begin with the Hamiltonian. Assuming only near- 
est neighbor interaction, it is given by H — Y^iLi 1 
where i's denote the sites in increasing order, and N is 
the total number of sites. H^j+i can be written as 
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FIG. 1: (Color online) Phase boundaries obtained by MF 
(black solid lines) and DMRG (red dotted lines). Upper- left: 
spin-spin correlation for the (0, 1,0) point. Upper-right: spin- 
spin correlation for 23 Na. (N = 60). The location of 23 Na, 
83 Rb, and 87 Rb within this parametrization are indicated by 
blue square, triangle, and circle respectively. 



where Ps,i,i+i denotes projection operators for sites i and 
i + 1 onto a state with total spin 5". Within second or- 
der perturbation theory in the hopping t between nearest 
neighbors, 65 = —At 2 /Us where Us is the Hubbard re- 
pulsion for two particles with spin S on the same site. 
Us is proportional to as, the s-wave scattering length 
in the spin S channel. For one-particle per site to be 
stable, we then need Us > for all S = 0,2,4, and 
hence €0,62,64 < 0. Within DMRG, however, it is more 
convenient to explicitly express -ff^i+i in terms of spin-2 
operators S,;, resulting 
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where ai = -l/3e -20/21e2+l/48e 4 , a 2 = -17/180e - 
l/9e 2 + l/40e 4 , a 3 = +l/45e - l/18e 2 + l/180e 4 , and 
a 4 = +l/180e - l/126e 2 + l/2520e 4 Qg^. Since it is 
crucial to identify the total spin for g round and excited 
states, we use non-Abelian DMRG [21| which allows us to 
find the lowest energy state in different total spin sectors: 
Stot = 0, 1, 2, 3, 4, and N. In this work we take N up to 
60 with open boundary condition (OBC), the number of 
states kept m = 400 and perform more than 20 sweeps to 
ensure the convergence. We note that though Eq. ([2} is 
more useful for numerical computation, it is more conve- 
nient to discuss the physics using Eq. ^ , which we shall 
do below. 

Since the phase diagram cannot depend on the overall 
energy scale, we plot the phase diagram in terms of the 
new variables (x , x 2 , £4) = (eo, e 2 , e 4 )/(eo + e 2 + 64). By 
definition < 2:0,2,4 < 1 and xq + x 2 + x± — 1. It is there- 
fore convenient to present our results using ternary phase 
diagram, where xo,2,4 are depicted as positions in an equi- 
lateral triangle. We place the points (0,0,1), (1,0,0), 



(0, 1,0) at the lower left, lower right, and upper corner 
of the triangle. The coordinates (xq, x 2l £4) at a given 
point within the triangle are to be read-off by drawing 
parallel lines towards the edges and then reading off the 
intersections. Our main result is as shown in Fig. [1] 

To describe and understand the results, we begin by 
giving the mean-field phase diagram as a useful reference. 
As mentioned, this can be found in direct analogy with 
the Bose-condensed case [ll|, 12 1 (see also [l3, 14, 17 1). 
The stability region for each phases are (in the form of 
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These are plotted as black-solid lines in Fig. [TJ F, P, and 
C occupy respectively the lower left, lower right, and up- 
per part of the triangle. Since the mean-field energies of 
the phases must be linear functionals of es's, the phase 
boundary are straight-lines. They all originate from the 
point where xq = x 2 — 24 = 1/3, the center of the 
triangle, where all three phases are degenerate. The in- 
tersection points at the edges are: FP: (17/24,0,7/24), 
PC: (10/17,7/17,0), CF: (0,1/2,1/2). 

Now we proceed to describe our phase diagram based 
on DMRG. Before we discuss each phase in detail, we 
summarize our finding as follows: The polar phase is re- 
placed by the dimcrized phase, while the cyclic phase is 
replaced by the trimerized phase. The FD and DT lines 
(no longer straight-lines) arc shifted away from the FP 
and PC lines but the FT line remains indistinguishable 
with the CF line. The phase diagram is easiest to un- 
derstand for the lower left region where X4 > xo. 2 . It 
is clear that the system would like to acquire the largest 
possible total spin. For N sites, the total spin is S = 2N. 
For this state, any two neighboring sites have total spin 
4 and the bond energy is 64, the smallest possible value. 
The ground state is thus 4iV + 1 fold degenerate. From 
our DMRG calculations wc find that the stability region 
for this state F actually extends slightly beyond X4 > x . 

For the rest of the phase diagram, it is more convenient 
to first consider finite N. Later, we shall mention how 
these pictures are modified as N — > 00. 

We find that the lower right region of the triangle is 
occupied by the dimcrized phase, where the unit cell is 
doubled, and the system is in a non-degenerate singlet 
state. In upper-right inset of Fig. Q] we provide an exam- 
ple for the spin-spin correlation function between neigh- 
boring sites, where the parameters of 23 Na are used. The 
correlation function clearly shows a "strong-weak" dimcr 
pattern. This phase is in direct analogy with the spin- 
1 case [22j]. The ground state for N ~ 2p sites can be 
most simply understood by first considering the case of 
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FIG. 2: (Color online) Excitation gaps to first few excited 
states with total spin S for 23 Na (upper plot) and (0, 1,0) 
point (lower plot) when a dimerization (trimerization) super- 
lattice of strength A2(As) is imposed. Open symbols: N = 60; 
Full symbols: lowest excited state for N — > do. 



xq 1 and imposing an artificial dimerization super- 
lattice where alternate bonds are weakened by a factor 
< A 2 < 1 (i.e., H 2l .2i+i -> A 2j ff 2 i,2m)- When A 2 = 0, 
the system breaks into N/2 subsystems, each consist- 
ing of only two interacting sites. For xq ~ 1, these two 
sites form a singlet with a finite gap to the first excited 
state(s), and the system is maximally dimcrized. As one 
increases A 2 , one expects that the gap should decrease 
gradually. In Fig. (left panel) we plot the gap for 
the first two excited states as a function of < A 2 < 1. 
We observe that the gap decreases monotonically with 
increasing A 2 but never vanishes, indicating the original 
ground state at A 2 = 1 is adiabatically connected with 
that at A 2 = and hence indeed dimcrized. 

Next we focus on the upper corner of the triangle, 
where we find the trimerized phase. For N = 3p, we 
find that the ground state is a spin singlet, with a fi- 
nite gap to the first excited state. In upper-left inset of 
Fig. [T] we provide an example for the spin-spin correla- 
tion function between neighboring sites for the case of 
X2 = 1- It shows a "strong-strong- weak" pattern, indi- 
cating the ground state is trimerized. To get a physical 
picture of this trimerized phase, it is helpful to consider 
a three-spins system in the limit where a; 2 ~ 1. For two 
interacting spins, the ground state is five-fold degener- 
ate, belonging to S = 2. For three spins (say 1, 2 and 
3) however, the ground state is a unique singlet: the 
total spin of any two spins can be 0,1,2,3,4, and only 
the spin-2 combination can be added to the third spin 
to form a singlet. By the same argument, this state is 
an eigenstate for both the operators H\2 and -ff 2 3, with 
eigenvalues e 2 . Hence the system has total energy 2e 2 , 
the lowest possible value. For system with TV = 3p sites, 
it is again helpful to impose a trimerization superlatticc 



where one every three bonds are weakened by a A3 factor 
(H 3l . 3l+1 -> A 3 i?3i,3i+i, < A 3 < 1). For A 3 -> 0, the 
system break up into iV/3 subsystems, each one with a 
singlet ground state just described and the whole system 
is maximally trimerized. The system at A3 = is hence 
gapped, and is stable towards increasing A3 from 0. In 
Fig. [2j3 (right panel) we show the gap as a function of 
A3. We find that the gap decreases monotonically with 
A3 but always remains finite, indicating that the ground 
state is adiabatically connected to the trimerized ground 
state of the A3 = limit. 

It is also instructive to study the behavior of the gap 
when a wrong superlattice is imposed. Consider a state 
in the dimerized (trimerized) phase. If one imposes a 
trimerization (dimerization) superlattice with strength 
As(A 2 ), the system will be converted into trimerized 
(dimcrized) phase as As(A 2 ) — > 0. As a result one ex- 
pects that there is a qualitative change in the ground 
state when As(A 2 ) is varied from 1 to 0, and the gap 
must vanish at a transition point. In Fig. [2] we also plot 
the gap as the function of A 2j 3 i when such a wrong super- 
lattice is imposed. We observe that the gap does vanish 
before A 2j 3 reaches zero, which provides an alternative 
confirmation for the nature of the phases. 

Now we consider N — > 00. For the dimerized phase, 
our description above remains valid, except that the mag- 
nitude of the gap becomes smaller. (Fig [21 full symbols). 
That is, the system remains gapped, unless a sufficiently 
strong "wrong" superlattice is imposed. For the trimer- 
ized phase, however, the situation is slightly different. 
Without any superlattice, the gap vanishes as 1/N as 
N ^ 00. That is, the system becomes gapless. Also, 
trimer order parameter defined by the appropriate sums 
and differences of the neighboring spin-spin correlation 
vanishes as N — > 00 (not shown). [23|. Nevertheless, 
for any finite A3, we find that the gap remains finite as 
TV — > 00. The gapped region as a function of A 2 .3 and 
1/N, as well as the gap dependence, are shown in Fig [3] 

It is natural to ask what are the expected phases of 
some of the available spin-2 elements. Using scattering 
lengths available in the literature 24|, we find that 83 Rb 
and 23 Na should be in the ferromagnetic and dimerized 
state respectively. For 87 Rb, however, the spin-spin cor- 
relation only shows a very weak trimcr-like pattern (not 
shown). We hence resort to use the method of imposing 
supcrlattices. We impose both trimerization and dimer- 
ization supcrlatticcs and calculate the gap as a function 
of A 2 ,3. We find that the gap monotonically increases un- 
der trimerization superlattice but decreases to zero under 
dimerization superlattice. This strongly suggests that 
the 83 Rb is indeed in trimerized phase. 

Spin-2 bosons in ID insulating lattice has been stud- 
ied theoretically before in Ref. [25(. Qualitatively, their 
picture agrees with ours for the lower left and lower right 
corners of the triangle, where they find ferromagnetic and 
dimcrized phase respectively. However, our results are 
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FIG. 3: (Color online) "Phase-diagram" (top panels) and en- 
ergy gaps (lower panels) as a function of N and A2,3. The 
red dotted lines indicate where the system becomes gapless. 
Left column, parameters according to 23 Na; right column, the 
point (0, 1,0). 



the corresponding system with A3 ; 2 = and then gradu- 
ally increasing it until it reaches 1. 

In summary we have studied the phase diagram of spin- 
2 bosons in ID optical lattice using non-Abelian DMRG. 
We identify three possible phases, namely ferromagnetic, 
dimerized, and trimerized, where the trimerized phase 
was not proposed in the literature. We demonstrate that 
by imposing proper superlatticcs and observing the exci- 
tation gap, the phase of interest can be identified or pre- 
pared based on the adiabatic connection of the ground 
state. Such a procedure can be implemented in cold- 
atom experiments and provide a simple scheme to test 
the phase diagram experimentally. 

We acknowledge inspiring conversations with M. A. 
Cazalilla. This research was supported by the National 
Science Council of Taiwan. 



qualitatively different for the rest of the phase diagram. 
Near the center of the triangle but excluding the ferro- 
magnetic regime, they claim to have nematic and cyclic 
phases. We do not find these states to be stable. Our 
result is also different in the region in the upper cor- 
ner. There, they proposed in fact two phases for this re- 
gion, which they named "cyclic" para-dimers and "para- 
dimerized" para-dimers. In their picture, every two sites 
form an effective spin 2 "para-dimers" , and these para- 
dimers either form a cyclic state, or more complicated 
composite objects which eventually have a cyclic order 
parameter. Trimcrization was not discussed in their pa- 
per. 

Lastly we discuss how our picture of the dimerized 
and trimerized state can be tested experimentally. There 
have now been many discussions on how to detect spin 
ordering for atomic gases in optical lattices [26|, [13 ■ We 
here, however, would like to point out a simple scheme 
particularly suitable for our states based on the adia- 
batic deformation of Hamiltonian by imposing supcrlat- 
tices as just discussed. Superlattices can easily be created 
[28j . and actually has already been utilized for prepara- 
tion and detection of certain spin-ordering states pjj . 
In our scheme, by applying an optical potential with 
two or three lattice spacing and changing the laser in- 
tensity for this optical potential, one can tune the pa- 
rameter A2,3 defined above since the tunneling ampli- 
tude t decreases with increasing potential barrier. A 
trimerizcd(dimerized) state is then characterized by a 
monotonic increase of the gap as a function of As(A2) 
and the gap closing before A2(Aa) reaches zero. Fur- 
thermore, by reversing the above argument, a trimeriza- 
tion(dimerization) supcrlattice can also allow us to pre- 
pare the trimerized(dimerizcd) state by first preparing 
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